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The ground state of a two-dimensional electron gas at Lan- 
dau level filling factors near v = 1 is a Skyrme crystal with 
long range order in the positions and orientations of the topo- 
logically and electrically charged elementary excitations of the 
v = 1 ferromagnetic ground state. The lowest energy Skyrme 
crystal is a square lattice with opposing postures for topo- 
logical excitations on opposite sublattices. The filling factor 
dependence of the electron spin-polarization, calculated for 
the square lattice Skyrme crystal, is in excellent agreement 
with recent experiments. 

PACS number 73.20.Dx 

The incompressible [0 ground state of a two- 
dimensional electron gas (2DEG) at Landau level filling 
factor v = 1 is ferromagnetic. (Here v = N/N^, is the ra- 
tio of the number of electrons to the orbital degeneracy 
of a Landau level; = AB/$ = A/{2i:f) where A is 
the area of the system, $o is the magnetic flux quantum, 
and B is the magnetic field strength.) Recently it has 
been shown that, in the limit of weak Zeeman cou- 
pling, the lowest energy charged excitations of this state 
are spin- textures known as skyrmions. fl . Skyrmions are 
the lowest energy topologically charged spin-texture ex- 
citations of the SU(2) non-linear a (NLcr) model which 
describes the long wavelength properties of isotropic fer- 
romagnets. The equivalence of physical charge and topo- 
logical charge in the present system is a consequence pg] 
of the quantum Hall effect jjj and is responsible for the 
dominating role of skyrmions in determining many phys- 
ical properties. 

Pioneering studies of skyrmions in the quantum Hall 
regime relied on NLcr models generalized (||| to account 
for Zeeman coupling to the spin and for the Hartree (elec- 
trostatic) interactions present when the charge density is 
non-uniform. States in the NLcr model are specified by 
a space-dependent unit vector which describes the local 
orientation of the spin magnetic moment. A skyrmion 
is characterized by the sign of its topological charge, by 
its size, and by the global orientation of the spin; with 
no Zeeman or Hartree coupling the energy of an isolated 
skyrmion is independent of all three. When Zeeman cou- 
pling is included the spin moment outside a skyrmion 
aligns with the field [m z = 1), the spin moment at the 
center of a skyrmion is oriented in opposition to the field 
(m z = — 1). The perpendicular component (m±) of the 
spin- moment, which must be non-zero as m z changes 



from —1 to 1, has a vortex at the skyrmion center. The 
global azimuthal orientation of m± is, importantly for the 
work described here, still arbitrary. Because the spin- 
moment is reversed in the interior of a skyrmion, Zee- 
man coupling favors small skyrmions. On the other hand 
the Hartree self- interaction energy of the skyrmion favors 
large skyrmions so that an optimal skyrmion size is estab- 
lished. H For Zeeman coupling strengths typical of phys- 
ical systems the estimated skyrmion size is comparable to 
microscopic lengths, motivating a microscopic approach. 
Recently |?]], using a microscopic Hartree- Fock approxi- 
mation, we confirmed the main predictions of the field- 
theory approach and obtained an estimate of the number 
of reversed spins for each (skyrmion) charge added to the 
system at v = 1. The quantitative agreement between 
this estimate and subsequent Knight shift measurements 
pj appears to provide incontrovertible evidence in favor 
of the exotic topological nature of the low-energy charge 
carriers in this system. In this Letter we report on a 
Hartree-Fock approximation study of the ground state of 
the 2DEG for v near 1, when interactions between the 
charged spin-textures, which we will refer to as skyrmions 
even when distorted by Zeeman and Hartree couplings, 
become important. 

The most important result of our calculation is shown 
in Fig. [l| where we compare our theoretical results for 
the spin-polarization P of the two-dimensional electron 
gas as a function of filling factor with experiment, j^] 
(P is the spatial average of m z .) Without interactions 
the spins are completely aligned by the Zeeman field for 
v < 1. The polarization is gradually reduced for v > 1, 
reaching zero for v — 2, since the Pauli exclusion prin- 
ciple forces the occupation of reversed spins. The large 
slope of P on either side of v — 1 is in sharp contrast with 
the non-interacting electron result and can be related to 
the number of reversed spins in an isolated skyrmion. 
HUJl]] Our calculations at finite \v — 1| were performed 
assuming that the ground state of the system is a crystal 
of charged skyrmion spin-textures. These states are anal- 
ogous to the Skyrme crystal states which arise || in stud- 
ies of dense nuclear matter using Skyrme's topological ex- 
citation model. The Hartree-Fock equations have many 
different solutions; all have long-range order in both the 
spatial arrangement of skyrmion centers and in the az- 
imuthal orientation of mi near each skyrmion center. 
Results are shown for two cases: i) a triangular lattice 
of skyrmions with identical orientations of m± near each 



site (TLF state) and, ii) a square lattice with oppositely 
directed orientations of m± on the two square sublat- 
tices (SLA state). (The initials A and F in the labels 
of these states suggest the obvious analogies with ferro- 
magnetic and antiferromagnetic states of spin systems on 
two-dimensional square and triangular lattices.) We see 
in Fig. ^ that the slope of the spin-polarization curve de- 
creases in magnitude as we move away from v = 1 in 
both cases, indicating that the skyrmions shrink as they 
become more dense. The shrinking is evidently much 
more rapid for the TLF state than for the SLA state. 
The results for P for the latter case are in quite rea- 
sonable agreement with experiment and indeed, as we 
discuss below, we find that the square Skyrme lattice has 
lower energy than the triangular Skyrme lattice except 
at very small \v — 1|. Comparison between theory and 
experimental results for P suggests that the ground state 
is a square Skyrme crystal until \v — 1| w 0.2. 

Our calculations are greatly simplified by working en- 
tirely within the lowest orbital Landau level and the 
Greens function equation of motion approach we have 
developed previously is readily adapted to the present 
situation. |nj For small enough | v— 1 1 , the skyrmions can 
be considered to be point particles and it is known from 
Madelung energy calculations for 2D lattices [jlljthat 
they will form a triangular lattice. Nevertheless, we find 
that by the time v = 1.1 the SLA state is lower than 
the TLF in energy by ss 30 times the amount by which 
the Madelung energy would favor the opposite energetic 
ordering of the states. The self-consistent spin-textures 
for the TLF and SLA states at v — 1.1 are illustrated in 
Fig. |^. The skyrmions are evidently more compact for 
the higher energy TLF state and this is the origin of the 
differences in spin-polarization between these states. 

To understand the differences between SLA and TLF 
states in greater depth it is useful to compare the con- 
tributions to the energy per electron from exchange, 
Hartree and Zeeman energies as shown in Fig. | for 
g/J-sB — 0.0l5e 2 /£. We express our results in terms 
of the energy (e qp ) increase per quasiparticle, N qp — 
\N — Nj,\, when the quasiparticles are created at fixed 
N by varying N^: 

_ u[e + (e 2 /l)^ + g^ B B/2] 
€qp - |1-„| W 

where e = E/N is the energy per electron, g is the g- 
factor for the 2DEG system, and {1b is the electron Bohr 
magneton. (The Hartree, exchange, and Zeeman contri- 
butions to e qp are given by terms of corresponding origin 
on the right-hand-side of Eq. [I];— (e 2 /£)^Tr/& is the ex- 
change energy per-particle at v — 1.) This convention 
for the energy per charge is implicit in the field theoreti- 
cal expressions and we choose it here because it will allow 
our Hartree- Fock results and results obtained in the field- 
theory description to be compared directly. In the limit 
of vanishing skyrmion density e qp is the so-called neu- 
tral quasiparticle energy. Il2| It follows from particle-hole 



symmetry that with this defintion antiskyrmion quasi- 
particle energies are smaller than skyrmion quasiparticle 
energies by g/isB. 

We first discuss the constituents of the energy for the 
Wigner crystal (WC) state which occurs near v = 1 in 
the limit of very strong Zeeman coupling. In this state 
a triangular (Wigner) crystal |lj| of single-reversed-spin 
quasiparticles is formed on the background of a full Lan- 
dau level of aligned spins. The reversed spins are local- 
ized to the greatest extent possible with the lowest Lan- 
dau level. For the WC state, the exchange and Hartree 
contributions to the energy per quasiparticle, e* p and e^L, 
can be calculated analytically jlj] in the limit of small 
"«>; 4 -> (e 2 /^)[(V8) 1/2 - U/16) 1/2 ] « 0.1835(e 2 /^) 
and eg (e 2 /£)( 7 r/16) 1 / 2 w 0.4431(e 2 /^), the Hartree 
self-interaction energy of a lowest Landau level quan- 
tized cyclotron orbit. The exchange energy, which is 
minimized when all spins are parallel, is nearly constant 
with increasing v qp until the quasiparticles begin to over- 
lap significantly while the Hartree energy is reduced at 
larger v qp due to spatial correlations of the quasiparti- 
cles. When the size of the quasiparticles (~ I for the WC 
state) is much smaller than the distance between quasi- 
particles the Hartree energy is reduced by the Madelung 
energy per particle [ [ill ; when the quasiparticles form a 

triangular lattice e rna d ~ —Q-7&21v qp 2 (e 2 /I). In Fig. || we 
have plotted the difference between the Hartree energy 
and the Madelung energy vs. v qp . For the WC state the 
remaining variation is due to the fact that the electrons 
are localized in cyclotron orbits and not at points. The 
Zeeman energy contribution to e qp is g^BB(S qp ) where 
(S qp ) is the number of reversed spins per quasiparticle. 
((S qp ) is smaller for antiskyrmions than for skyrmions by 
1.) For the WC state (S qp ) = 1 for skyrmions and for 
antiskyrmions. 

In the absence of Zeeman coupling a skyrmion quasi- 
particle has lower energy than a reversed spin quasipar- 
ticle because its charge is spread out over an arbitrarily 
large distance, reducing e^ p to zero while at the same 
time keeping spins locally nearly parallel so that 
increases only |,| to (e 2 /f)( 7 r/32) 1 / 2 w 0.3133(e 2 /^). 
When a finite Zeeman coupling is present, the skyrmion 
quasiparticle must shrink to reduce its Zeeman energy 
and the Hartree energy increases. For an isolated 
skyrmion (v qp — > 0) and g^sB = 0.015e 2 /^ we find us- 
ing the approach ||] of Ref. § that ef p = 0.2499(e 2 /f), 
= 0.2957(e 2 /f) and that (S qp ) = 3.83. In Fig. | we see 
that for the TLF state the Hartree energy quickly rises 
from the isolated skyrmion value toward the WC value 
and (Sqp) quickly decreases from the isolated skyrmion 
value to 1 as v qp increases. Evidently there is a strong 
short-range repulsive interaction between skyrmions in 
the TLF Skyrme state which quickly shrinks them into 
single-reversed-spin quasiparticles. The origin of this re- 
pulsion is an increase in exchange energy in the region 
between neighboring skyrmions with identical azimuthal 
orientations of m±. As we see in Fig. the azimuthal 



orientation of m± must then change by tt relatively 
abruptly along the line connecting neighboring skyrmion 
centers; in the field-theoretical description |||| the asso- 
ciated increase in exchange energy is described by a term 
oc |VtojJ 2 and the resulting repulsive interaction will be 
roughly proportional to m\ evaluated at the mid-point 
between skyrmion centers. Clearly this repulsive interac- 
tion is minimized if neighboring skyrmions have opposite 
azimuthal orientations for m± as illustrated for the SLA 
state in Fig. |[ (It seems clear that the lowest energy 
triangular Skyrme lattice will be a three-sublattice state 
with 27r/3 relative orientation but this state will, we be- 
lieve, be higher in energy than the SLA state unless v is 
very close to 1.) The exchange energy thus favors bipar- 
tite lattices for a Skyrme crystal. Our numerical results 
for the SLA state show that the skyrmions shrink rel- 
atively slowly as v qv increases; for this arrangement the 
Skyrme crystal is able to tolerate much tighter packing of 
the skyrmions and the rapid decrease in spin-polarization 
persists out to quite large values of v qv as shown in Fig. |l| 
and observed experimentally. Q 

The size of a quasiparticle in a Skyrme crystal is lim- 
ited by the lattice constant of the crystal until, with de- 
creasing v qp the lattice constant becomes substantially 
larger than the size of an isolated skyrmion quasiparti- 
cle. For weaker Zeeman couplings, isolated skyrmions 
are larger [^]|7|] and spin-polarization will decrease more 
rapidly for v very close to v — 1 ; for larger values of | v— 1 
where the skyrmion size is limited by inter-skyrmion in- 
teractions, the spin-polarization is much less dependent 
on the Zeeman coupling. For example we find for the 
TLF state at gii B B/(e 2 /£) = 0.0015 that P = 0.7 at 
v = 1.03, substantially smaller than the values reported 
in Fig. [I], while for v = 1.15, P = 0.7 again, nearly iden- 
tical to Fig. 0. 

From our calculations we can conclude that for typical 
qHbB values, the ground state of the 2DEG near v = 1 is 
composed of skyrmionic charged excitations of the ferro- 
magnetic v = 1 background. When these quasiparticles 
are sufficiently dilute they will crystallize because of the 
absence of kinetic energy. Although we have not com- 
pleted an exhaustive survey of possible Skyrme lattices, 
we believe that SLA Skyrme crystal will be lowest in en- 
ergy once the lattice constant is small enough for the gra- 
dient energy to overwhelm the relatively modest struc- 
tural preferences of the Hartree energy. It is difficult for 
us to estimate the range of filling factors over which the 
Skyrme crystal is favored over fluid states of the 2DEG. 
However, it seems unlikely that the energy competitions 
which determines skyrmion size will be very sensitive to 
the existence of long-range-order. We anticipate the ex- 
istence of a Skyrme liquid state, which will be lower in 
energy than maximally spin polarized liquid states and 
which will have a spin polarization quite close to that of 
the SLA state. The Skyrme liquid should have a typical 
coordination number closer to four than to six because 
of the dependence of the exchange energy on the relative 



azimuthal orientations of m± for nearby skyrmions. 
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tion between single-skyrmion states with adjacent values 
of (S z ). For typical systems we do not expect the dis- 
creteness of the allowed values for {S z } to be observable 
above ~ IK. In the Skyrme crystal skyrmion-skyrmion 
interactions will cause (S z ) to fluctuate on each site even 
at T = 0. 

FIG. 1. Variation of the spin polarization, P, as a func- 
tions of v, for SLA and TLF Skyrme crystal states and for 
different values of g* = g^BB/{e /£). The open circles are 
experimental data obtained for g* near 0.016, while the filled 
circles show data obtained for g* near 0.021. The spin polar- 
ization in the WC state is the same as that for independent 
electrons. 

FIG. 2. Two-dimensional vector representation of the x—y 
components of the spin density of a crystal of skyrmions as 
obtained in the HFA. The length of the arrows is proportional 
to the local magnitude of m± and their direction indicates the 
local orientation. The magnitude of m± is small for m z near 
-1, near the center of a skyrmion, and for m z near +1, outside 
of a skyrmion. (a) TLF state at v = 1.1 and g* = 0.015; (b) 
SLA state at v = 1.1 and g* = 0.015. 

FIG. 3. Zeeman (a) , Hartree (b) and exchange (c) contri- 
butions to the energy per quasiparticle for the Wigner crys- 
tal state (*) and for the TLF (triangles) and SLA (squares) 
Skyrme crystal states. The results for v = 1 were obtained 
analytically for the WC state and from isolated skyrmion cal- 
culations for the Skyrme crystal states. The Zeeman energy 
is expressed in terms of the average number of reversed spins 
per quasiparticle and the Madelung contribution has been 
subtracted from the Hartree contribution, as discussed in the 
text. 



